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The magnetic flux dynamics of type-II superconductors within the critical state regime is posed 
in a generalized framework, by using a variational theory supported by well established physical 
principles. The equivalence between the variational statement and more conventional treatments, 
based on the solution of the differential Maxwell equations together with appropriate conductivity 
laws is shown. Advantages of the variational method are emphasized, focusing on its numerical 
performance, that allows to explore new physical scenarios. In particular, we present the extension 
of the so-called double critical state model to three dimensional configurations in which only flux 
transport (T-states), cutting (C-states) or both mechanisms (CT-states) occur. The theory is ap- 
plied to several problems. First, we show the features of the transition from T to CT states. Second, 
we give a generalized expression for the flux cutting threshold in 3-D and show its relevance in the 
slab geometry. In addition, several models that allow to treat flux depinning and cutting mech- 
anisms are compared. Finally, the longitudinal transport problem (current is applied parallel to 
the external magnetic field) is analyzed both under T and CT conditions. The complex interaction 
between shielding and transport is solved. 

PACS numbers: 74.25.Sv, 74.25.Ha, 41.20.Gz, 02.30.Xx 



I. INTRODUCTION 

The investigation of the macroscopic magnetic prop- 
erties of type-II superconductors in the mixed state is 
already a classical subject. The essential physics behind 
the collected vast phenomenology has been well known 
for decades fii^ and may be basically analyzed in terms 
of interactions between the flux lines themselves (lattice 
elasticity, line cutting), and interactions with the under- 
lying crystal structure (flux pinning). 

For many purposes, it happens that the mesoscopic de- 
scription may be further simplifled by means of appropri- 
ate material laws for the coarse-grained fields, i.e.: mag- 
netic induction B = (b), electric current density J = (j) 
and electric field E = (e). Averages are supposed to be 
taken over a volume containing a big enough number of 
vortices. Along this lines, a brilliant proposal, the so- 
called critical state (CS) theory, originally introduced by 
C. P. Bean^ has been widespread used. Such a model 
allows to capture the main features of the magnetic re- 
sponse of superconductors with pinning at low frequen- 
cies and temperatures, through the minimal mathemat- 
ical complication. In its simplest form, the CS theory 
involves to solve Ampere's law dB/dx = hqJ with some 
prescription for the current density (J = ± or 0), and 
under continuity boundary conditions that incorporate 
the influence of the sources. Being a quasi-stationary 
approach, the CS is customarily stated without an ex- 
plicit role for the transient electric field and the remain- 
ing Maxwell equations. Nevertheless, the recognized pre- 
diction power of the theory is not accidental. Although 
veiled, the role of E and Faraday's law is of great im- 
portance for its soundness. Thus, we recall the above 
statement may be related to an almost vertical E{J) law. 



idealized by the graph E' = for J < Jc and E' ^ oo for 
J > Jc- Here, E stands for the induced electric field ow- 
ing to variations of the fiux density, and Faraday's equa- 
tion (in fact Lenz's law) is implicitly used by selecting 
the actual value ± Jc or that minimizes flux variations, 
when solving dB/dx = (IqJ along the process. 

In the case of ideally one-dimensional problems, i.e.: 
infinite cylinders and slabs in parallel field configuration, 
the previous statements lead to the prediction of the ob- 
servable physics without ambiguity. As recalled, one is 
making a correct use of the Maxwell equations for a con- 
ducting system that experiments a sharp transition in 
terms of the current density. Physically, the material 
law relates to the vortex pinning phenomenon, that is 
macroscopically described by the average pinning force 
constraint |J x B| < Fp. Straightforwardly, this gives 
place to the depinning threshold limitation J± < Jc± (B) . 
In addition, one has to consider a high flux flow dissipa- 
tion when the limitation is exceeded (as related to the 
mathematical condition E — > oo). We call the readers 
attention that for the simplified geometries, the equality 
J = Jx (J is locally perpendicular to B) is automatically 
fulfilled. 

The extension of the above ideas to non-idealized sys- 
tems is not a closed subject yet, it is of utter impor- 
tance for the understanding of the experimental facts, 
and constitutes the main motivation of this article. In 
brief, and following the spirit of Bean's model, we con- 
sider the question of identifying a theory that allows a 
general description of the low frequency electrodynamics 
of hard superconductors with the least conceptual and 
mathematical load. 

From the physical point of view, the allowance of non- 
parallel flux lines leads to consider a new threshold for 
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the current density, now related to the disorientation of 
adjacent vortices. This was aheady remarked by J. R. 
Clem and A. Perez-Gonzalea^ who analyzed the relation 
J|| = Hda/dx with a the angle characterizing the flux 
lines orientation at a given depth within a superconduct- 
ing slab. Thus, stemming from the fact that a max- 
imum angle gradient is allowed between vortices so as 
to avoid cutting and recombination, one has to consider 
•^11 < Jc\\{B):^ Remarkably, the above mentioned au- 
thors showed that both the flux depinning and cutting ef- 
fects may be treated in a generalized CS framework. The 
upgraded theory (double critical state model or DCSM) 
has been applied with high success since the early 80 's 
for the understanding of many experiments in which 
both J|[ and Jj_ are involved. In addition to the men- 
tioned static threshold conditions, the DCSM general- 
izes the one-dimensional E{J) graph for quasi-stationary 
processes in terms of the natural concepts £'|[(J|[) and 

From the mathematical point of view, the DCSM has 
been mainly applied to experiments with rotating mag- 
netic field components, still parallel to the surface of 
large samples, i.e.: two-component fields are allowed, but 
only one independent variable is considered. However, 
some aspects that have appeared linked to the investi- 
gation of new materials (small crystallites and high-Tc 
films), as well as the refinement of previous studies re- 
quire more specialized statements, so as to include fi- 
nite size effects, sample inhomogencity, anisotropy, etc. 
Along this line, some recent advances have to be quoted. 
On the one side, it has been shown that finite size ef- 
fects for thin samples may be treated by composition of 
quasi one-dimensional statements)^ In principle, this idea 
would allow to include both the Jci_ and Jc\\ limitations 
but, to the moment, it has been exploited in the limit 
Jc± <C Jc\\ ■ On the other hand, truly two-dimensional 
configurations in which some symmetry property allows 
to assume J _L B have also been solved by numerical 
methodsi^i^ More recently, as a remarkable advance to 
be mentioned, the mathematical structure of the DCSM 
solution in the three dimensional case has been described 
and obtained variationally for some examplesJ^ 

In this work, new perspectives on the application of 
variational methods^i will be presented. Contrary to 
some recent claim about their restricted scopeji they will 
be shown to be equivalent to the more conventional dif- 
ferential equation statements for solving CS problems. 
Moreover, in many instances, our solution will be used 
to extend previous results toward unexplored physical 
scenarios. In particular, the influence of the parame- 
ter ratio X ^ Jc\\/Jc± will be quantified in three di- 
mensional systems. Also, allowed by the capability of 
the theory, we introduce a new critical angle gradient 
threshold, that generalizes the two dimensional concept 
da/dx = J\\/H < Kc- On the other hand, the so-called 
longitudinal transport problem, i.e.: a situation in which 
a magnetic field is applied along the direction of the 
transport current will be studied in a 3-D configuration. 



As a central physical result of our paper, it is shown that 
the variational method naturally distinguishes between 
the inductive and potential parts of the background elec- 
tric field in the CS problems. It will be shown that the 
incorporation of the physical idea of the direction of the 
electric field is straightforward by combining Gauss' law 
and the variational statement. 

The paper is organized as follows. In Sec. [TTl the phys- 
ical background of the general CS concept is introduced. 
The underlying approximations pi A[) . the validity of the 
associated material law (jll Bp . and the justification of a 
variational statement (jll C[) are described in detail. In 
Sec. mil we give a number of explicit examples related to 
the application of the general CS theory to three dimen- 
sional systems. Specifically, we consider various magnetic 
processes for an infinite slab with a penetrating magnetic 
field of the form (Hx, Hy, Hz), and described by differ- 
ent models hosted in our theory. It will be shown that 
our results fully coincide with alternative formulations, 
when comparison is allowed. Finally, Sec. IIVI is devoted 
to discuss the main results of this work. 



II. GENERAL CRITICAL STATES: THEORY 

This section is devoted to introduce the theoretical 
background that justifies the critical state concept as 
a valid constitutive law for superconducting materials. 
First, recalling that the CS must be considered an ap- 
proximation within the magnetoquasistationary (MQS) 
regime of the time-dependent Maxwell equations, we will 
discuss on the related limits for the whole set of elec- 
tromagnetic physical quantities. Then, we will present a 
thorough discussion about the representation of the CS 
as a J(E) law for a perfect conductor with restricted cur- 
rents, in the MQS regime. Finally, the variational state- 
ment of the CS is stated. We will show that the varia- 
tional principle introduced in previous workiiii^ is fully 
equivalent to the more usual formulation in terms of dif- 
ferential equation statements, in time-discretized form. 
Along this line, some recent claims about the limited 
scope of variational formulations^ have to be reconsid- 
ered. 



A. MQS approximation and its consequences 

Let us first concentrate on the physical implications re- 
lated to the MQS approximation within the critical state 
theory. Recall that, in general, the dynamical behavior 
of the macroscopic electromagnetic fields is determined 
by the Maxwell equations accompanied by material con- 
stitutive laws, H(B), D(E), and J(E). Thus, Faraday's 
and Ampere's laws represent a coupled system of time 
evolution field equations 

5tB = -V X E , 9tD V X H - J . (1) 
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FIG. 1: (Color online) Schematics of the time dependence of 
the electromagnetic fields within the MQS approximation. A 
ramp in the magnetic field is induced by the external excita- 
tion, within the interval + As a consequence of a very 
fast diffusion (elevated fiux flow resistivity), the electric field 
quickly adjusts to a constant value along the interval. When 
the magnetic field ramp stops, E goes back to zero again. 
The re-adjusting vertical bands are considered a second or- 
der effect and allow for charge separation and recombination, 
according to the specific E(J) model. 



Taking divergence in both sides of each, and recalling 
integrability (permutation of space and time derivatives) 
leads to the additional conditions 

dt{V-B) = Q , 9t(V-D) + V-J = 0. (2) 

Within this picture of the electromagnetic problem, the 
remainig Maxwell equations can be interpreted as "(spa- 
tial) initial conditions" for Eq.(I2]), that define the exis- 
tence of conserved electric charges, i.e.: 

V-B(t = 0) = , V •D(t = 0) = = 0) . (3) 

Eqs. ([1]), upon substitution of H, D and J through 
the constitutive laws, and with appropriate initial con- 
ditions, uniquely determine the evolution profiles B(r,t) 
and E(r, t). 

In this paper, as it is customary in hard superconduc- 
tivity, we will consider B = /ioH as a valid approxima- 
tion related to neglecting the equilibrium magnetization 
of the flux line lattice. On the other side, for slow and 
uniform sweep rates of the magnetic sources, the tran- 
sient variables E, D and p are small, and proportional to 
B, whereas B, E and p are negligible. Thus, the main 
hypothesis within the MQS regime is that the displace- 
ment current densities 9tD are much smaller than J in 
the bulk, and vanish in a first order treatment (see Fig. 
[T]). This causes a crucial change in the mathematical 
structure of the Maxwell equations: Ampere's law is no 
longer a time evolution equation, but becomes a purely 
spatial condition. It reads 

V X B ^ /^oJ , (4) 



with approximate integrability condition V ■ J ~ 0. In 
the MQS limit, Faraday's law is the unique time evo- 
lution equation. Then, making use of the conductivity 
law through its inverse function E(J), one can find the 
evolution profile B(r, t) from 

dtB = -V X E(^oJ =i V X B) . (5) 

We want to mention that the formulation in Eq. ^ is 
definitely the most extended one. However, the possibili- 
ties of formulationai^ and J— formulationsii in which 
the dependent variables are the fields E and J have also 
been exploited by several authors. Also, a vector po- 
tential oriented theory (A— formulation) has been issued 
recentljii^, that is a promising path for the investigation 
of 3-D problems. 

Let's point out two relevant consequences of the MQS 
limit. 

1. The constitutive law D(E), which is not used in 
Eq. (ini), plays no role in the evolution of the mag- 
netic variables B and J; the magnetic "sector" is 
uncoupled from the charge density profile because 
the coupling term (charge recombination) has dis- 
appeared. 

2. Only the inductive component of E (given by 
V X Eind — — B, V • Eind = 0) determines the evo- 
lution of B (Faraday's law). The conducting law in 
its inverse formulation E(J) presents some ambi- 
guity, as far as two different material laws related 
by E2(J) = El (J) + V'I'(J) determine the same 
magnetic and current density profiles. 

Going into some more detail, whereas for the complete 
Maxwell equations statement, the potential component of 
the electric field (V x Epot = 0, cqV • Epot ~ p), is cou- 
pled to B and Eind through the D term, within the MQS 
limit it is irrelevant for the magnetic quantities. In fact, 
one can include the presence of charge densities without 
contradiction with the condition V- J ~ by means of in- 
homogeneity or non-linearity in the E(J) relation. Then 
one has V-J = O^V-E = 0. The charge density p can 
be understood as a parameterized charge of static char- 
acter, as far as p is neglected. As indicated above, once 
integrated the magnetic variables, and computed V • E, 
one has the freedom to modify the electrostatic sector if 
necessary, by the rule E(J) + V^, while still maintaining 
the values of B and J. This invariancc can be of practical 
interest, as far as the "electrostatic" behavior in the CS is 
still under discussion, because of the inherent difficulties 
in the direct measurement of transient charge densities. 
Recent advances have to be quoted^^, but they are still 
based on the analysis of Ejnd and some ansatz on the 
direction of the electric field. To be specific, E || J is 
assumed in that work. 
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FIG. 2: (Color online) Schematic representation of the CS 
E{J) model. The electric field arises when some critical con- 
dition for the volume current density is reached (Jc in this 
1-D representation). Corresponding to the MQS approxima- 
tion in Fig. [1] the electric field instantaneously increases to 
a certain value, determined by the rate of variation of the 
magnetic field, and then goes back to zero. 



B. Material law: the critical state 

Now, we will be more explicit about the material law 
J(E) that dictates the magnetic response of a supercon- 
ducting sample in the critical state and for a given exter- 
nal excitation. For simplicity, we start with an overview 
of the material law for 1-D systems (infinite slabs or 
cylinders with the external field applied along symmetry 
axis). The physical concepts will be eventually general- 
ized to 3-D. 



1. 1-D Critical States 

For our purposes, it is sufficient to recall that the ba- 
sic structure of the CS relates to an experimental graph 
within the {V, /} plane that basically contains two re- 
gions: 

1. —Ic < I < Ic with perfect conducting behavior, 
i.e.: V = and dj = 0. 

2. For I ^ Ic, the curve is characterized by a high 
diV slope (and antisymmetric for / < —Ic)- Fur- 
ther steps, with / increasing above the critical value 
/c, i.e., the eventual transition to the normal state, 
may be neglected for slow sweep rates of the exter- 
nal sources, which produce moderate electric fields. 

Within the local description level, different models have 
been used for the corresponding E ^ J graph, the most 
popular being 

1. The power law model: E = a sgn(J) (| J|/ Jc)"", 
with a a constant and n high. 



2. The piecewise continuous linear approximation: 
£; = for I J| < Jc, and E = f3 sgn(J)(| J| - Jc) 
for I J| > Jc, (3 having a high value. 

This model and the previous one present a small 
dependence on the sweep rate, as far as different 
values of E give way to a slightly different J. 

3. Bean's model: constant J for i? = 0, and J = 
sgn{E)Jc for E ^0 (see Fig. ^. 

This is the simplest model, without sweep rate de- 
pendence because only the sign of E enters the the- 
ory 

Bean's model captures the main features of the CS and 
has been widely used with very good performance since 
its proposal in the early sixtiesi^ On the other hand, we 
notice that it may be obtained from the other represen- 
tations: it corresponds to the limiting cases n oo and 
f3 —>■ oo respectively, as the reader can easily check. Con- 
sidering the above ideas as the essential hypotheses of the 
CS theory, the 3-D critical state model will be formulated 
upon its generalization. The well known experimental 
evidence of a practical sweep rate independence for mag- 
netic moment measurements (unless for high frequency 
AC sources or at elevated temperatures) reinforces this 
simple model as a valid tool in the CS. 

In some treatments, the first or second models arc im- 
plemented, in order to transfer a full E(J) law to the 
Maxwell equations. On the other hand, being rate in- 
dependent. Bean's model is no longer a time-dependent 
problem, but a path-dependent one, i.e.: the trajectory 
of the external sources Hg uniquely determines the mag- 
netic evolution of the sampleJ^ This makes an important 
difference when one compares to more standard treat- 
ments, as far as Faraday's law is not completely deter- 
mined from the path. Strictly speaking, one has 

AB === -V X [EAt] , (6) 

with At (and therefore |E|) gauged by the evolution of 
the external sources. In other words, the absence of an 
intrinsic time constant gives way to the arbitrariness in 
the time scale of the problem. On the other hand, the 
magnitude |E| is not relevant for the B and J profiles. In 
fact, in the applications of Beans's model, Faraday's law 
is not strictly solved. It is just the sign rule (the vectorial 
part of the material law), that is used to integrate Am- 
pere's law. Notice that such sign rule corresponds to a 
maximal shielding response against magnetic vector vari- 
ations, and thus, determines the selection of J = ±Jc. 

Notice that, by symmetry, in the 1-D problems one has 
J II E and both quantities are orthogonal to B. Thus, 
at a basic level, the 1-D CS concept is grounded on the 
existence of pinning forces that act as a barrier against 
flux flow. The physical threshold related to a maximum 
value of the force balancing the magnctostatic term J x B 
gives place to the concept of maximum (critical) current 
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density, and thus to the law 

^ sgn{E^)J,^ for E^^O. (7) 

Here, E± stands for the component of E along the direc- 
tion B X (J X B). 

2. General (3-D) Critical States 

Let us now see how the above ideas may be trans- 
lated to a 3-D scenario from the fundamental point of 
view. The main issue is that, in general, the paral- 
lelism of E and J and their perpendicularity to B are 
no longer warranted. Then, a sign rule does not suffice 
for determining the solution. A vectorial rule is needed 
and attention must be paid to its mathematical consis- 
tency, as well as to the physical significance. In previous 
work,— li^ii^ we introduced a geometrical concept that 
may be of much help when discussing the idea of a gen- 
eral critical state theory. There must be a region Ar 
within the J-space (possibly oriented according to the 
local magnetic field B, and/or also depending on |B| and 
r) such that non-dissipativc current flow occurs when the 
condition J e Ar is verified (see Figs. |3l E]). This con- 
cept, together with a very high dissipation when J is 
driven outside Ar by some non vanishing electric field, 
suffice to determine the relation between the directions 
of J and E. This may be done according to the following 
argument. 

7k- In the critical state, the forces arising to avoid the 
flux flow (or whatever) dissipation mechanism act 
against the exit of J from a region Ar with a very 
high slope. However, the evolution within Ar is 
that of a perfect conductor. In the limit of an infi- 
nite barrier, the reaction is perpendicular to the 
boundary of Ar (denoted by dA^). Thus, start- 
ing from an initial configuration with J £ Ar, and 
under the action of a transient electric field, the 
vector J quickly touches and/or shifts along the 
boundary, until a point is reached where the con- 
dition E _L 9Ar is fulfilled. Owing to the perfect 
conductivity condition dtJ oc E no further evolu- 
tion can occur (see Fig. [3] and ReffTol) . Faraday's 
law will be eventually the key for determining the 
actual point. 

Recall that the above rule can be expressed as the con- 
dition of maximal J.E projection for J G Ar. Recall also 
that the fundamental property already discussed for 1- 
D systems is verified: still the modulus |E| is irrelevant. 
On the other hand, notice that the sign rule is nothing 
but the 1-D particular case of the general maximum pro- 
jection rule. The mathematical consistency is therefore 
satisfied. Let us now discuss on the physical soundness 
of the theory. This aspect is apparently related to the 
selection of the region Ar appropriate for the problem 
under consideration. 




FIG. 3: Top: schematic representation of the relative orienta- 
tions of the local magnetic field H and electric current density 
J. The current is decomposed into its parallel and perpendic- 
ular components, i.e.: J = Jy + J±. Also sketched is the di- 
rection of the magnetic field at some neighboring point, at an 
angle The vectors H, H' and J do not necessarily he at the 
same plane. Bottom: the perfect conducting region within the 
plane perpendicular to the local magnetic field. An induced 
electric field is shown. Initially (Jo), the high dissipation re- 
gion is touched, but almost instantaneously J shifts along the 
boundary, reaching a point where the condition E _L dAr is 
fulfilled. Anisotropy within the plane is allowed. 



The simplest assumption that translates the CS idea 
to 3-D situations was already issued by Bean in Ref. [20l . 
It has been called the isotropic CS model and generalizes 
1-D Bean's law to 

J = Jc E if S ^ , (8) 

i.e., the region Ar becomes a sphere. This model has 
been used by several authors; in spite of its mathematical 
simplicity, a remarkable predictive power for reproduc- 
ing a number of experiments with rotating and crossed 
magnetic fields has been noticed )^^'^^i^^'^^ at least qual- 
itatively. It lacks, however, a solid physical basis. In 
any case, one could argue that statistical averaging over 
a system of entangled flux lines within a random pinning 
structure might be responsible for the isotropization of 

Ar. 

As stated above, to the moment, the most general the- 
ory for CS problems, formulated in terms of a well ac- 
cepted physical basis is the so called double critical state 
modelA (DCSM). In brief, this theory assumes two differ- 
ent critical parameters, Jc\\ and (see Fig. [J]) acting 
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FIG. 4: Geometric interpretation of the Critical State behav- 
ior for the DCSM case. J is constrained to the boundary of 
a rectangular region. T, C and CT states are related to the 
horizontal and vertical sides, and to the corners. Two models 
in which the corners of the DCSM region have been smoothed 
are also shown. 



as the thresholds for the components of J parallel and 
perpendicular to B respectively. As stated above, Jci_ 
relates to the flux dcpinning threshold induced by the 
Lorentz force on flux tubes, while the additional J^w is 
imposed by a maximum gradient in the angle between ad- 
jacent vortices before mutual cutting and recombination 
occursi^ The DCSM may be expressed by the statement 
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(9) 



Within the DCSM, the region Ar is a cylinder with its 
axis parallel to B, and a rectangular longitudinal sec- 
tion in the plane defined by the unit vectors B, J_l (see 
Fig. 2]). The edges of the region Ar introduce a criterion 
for classifying the CS configurations into : (i) T-states 
where the flux dcpinning threshold has been reached (J 
belongs to the horizontal sides of the rectangle), (ii) C- 
states where the cutting threshold has been reached (J 
belongs to the vertical sides of the rectangle), and (iii) 
CT-states where both J|| and J± have reached their crit- 
ical values (corners of the rectangle). 

Notice that Jf,|| and Jcj. are determined from different 
physical phenomena, and their values may be very dif- 
ferent (in general Jc\\ > Jc± or even Jc\\ ^ Jc±)- Never- 
theless, the coupling of parallel and perpendicular effects 
is suggested by experiments^^ and, for instance, may be 
included in the theory by the condition Jc\\ ~ KBJc±_ 
with K a material dependent constant. Recalling that 
the mesoscopic parameters Jc are related to averages over 
the flux line lattice, interacting activation barriers for the 
mechanisms of flux dcpinning and cutting are expected 
and this may give place to deformations in the boundary 
dAr- Then, the theory should be able to host different 
regions as the ones depicted in Fig. [H or the situation 
suggested in Ref. (Fig. 1). Regarding that proposal. 



we want to emphasize that the statement qualifying our 
maximum projection rule (i.e.: max J-E <^4> E _L 9Ar) as 
physically incorrect in some cases, because the direction 
of E should be different, must be reconsidered. In fact, as 
exposed above, in the CS theory the material law allows 
some ambiguity described by an arbitrary scalar function 
in the form of an additional potential term, i.e.: the rule 
E = Ecs + V(f> docs not affect the magnetic variables. 
Here, Ecs represents the electric field obtained from the 
maximal shielding rule, and E a possible modification in 
order to adjust the electrostatic sector, f.i. the scalar 
condition E • B = (E ± B) for the T-states. Thus, 
the maximal shielding rule can be easily complemented 
with an additional equation V<I> • B = —Ecs ■ B for the 
new variable <f> that allows to link our theory to many 
other models, expressed by an E(J) law. Eventually, not 
only the magnetic sector but also the electrostatic one 
would coincide in both formulations. In physical terms, 
one can state the theory as follows. First, the magnetic 
response is described by the Critical State law, that de- 
termines the relative orientation of Ecs and J through 
the maximum projection rule. If required, the orienta- 
tion of the full electric field may be tuned by considering 
the electrostatic charges. 



C. The variational principle: a general treatment 

The CS theory has been formulated as a minimization 
principle by several groups (see f.i.: Ref s fllllT^l2^l25H26h . 
In our case^iiii^ an Optimal Control^ variational state- 
ment was introduced for dealing with general CS prob- 
lems. In previous work,— a number of experimental facts 
were discussed in terms of different choices for the control 
region Ar. Here, we want to emphasize that the varia- 
tional principle is not to be associated to any particular 
model, i.e., an arbitrary selection of the region Ar is al- 
lowed. 

Below, we summarize the main features of the formu- 
lation. Its full equivalence to the approach based on the 
differential equation statements will be shown. The hy- 
potheses about constitutive relations and the sweep rate 
independence made in the earlier analysis will be obvi- 
ously maintained. The principle is based on a discretiza- 
tion of the path followed by the external sources, that is, 
it is an approximation to the continuous evolution whose 
accuracy increases as the step diminishes. 

Let us consider a small path step, from some ini- 
tial profile of the magnetic field Bn(r) to a final pro- 
file Bn+i(r), (define AB = Bn+i — Bn, and also the 
corresponding Jn(r) and Jn-(-i(r)). Both configurations 
can be considered to be connected by a stationary pro- 
cess, i.e., we perform a small linear step AB, such that 
Bn+i = Bn + sAB, s G [0, 1]. The initial condition ful- 
fills Ampere's law V x Bn = /xnJn, as well as V • Bn = 0, 
V • Jn = 0. As shown in Ref. Il2l . maximal shielding can 
be implemented by imposing the minimization of the step 
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variation for the magnetic field profile integral, i.e. 



^[B„+i(-)] = Min 



(10) 



Recall that minimization must be performed under 
the restrictions on the final profile: (i) Ampere's law 
V X Bn+i = ^oJn+1, and (ii) Jn_|_i(r) e Ar. This is 
a minimization problem within the variational calculus 
framework (integral functionals of unknown fields and 
their derivatives) with constraints, that can be analyzed 
with the tools of the Optimal Control theoryi^ 

Following the usual Lagrange multipliers method, we 
build a "Lagrangian" 



L = J(AB)2 + p • (V X Bn+l - Jn+l) 



(11) 



that enforces Ampere's law. In fact, the Euler-Lagrangc 
equations become 

V X Bn+i - J„+i = (12) 

for arbitrary variations 5p of the multipliers, and 



V X p = -AB, 



(13) 



for arbitrary variations 6(Bn+i). 

The second condition identifies p with — AA (recall 
that V X AA ~ AB). Then, one gets the critical state 
electric field EcsAt = — AA = p. Concerning the "pa- 
rameter" Jn+i, as far as it is not allowed to take arbi- 
trary values, we cannot impose arbitrary variations as it 
is customary for the typical stationarity condition of the 
Euler-Lagrange equations. Instead, an Optimal Control- 
like Maximum principle must be usedf^ The minimum of 
the Lagrangian must be sought within the set of current 
density vectors fulfilling J G Ar, i.e.: Jn-i-i is determined 
by the condition 



Min{L} |jeA, 



Max {J • p} |jeA, . (14) 



Notice that the maximal shielding condition is equiva- 
lent to the maximum projection rule, i.e.: the orthog- 
onality condition of the electric field direction with the 
surface of Ar previously discussed (Sec lII B| ) is recovered. 
Notice also that Ampere's law is imposed [Eq. (HH)] 
through the Lagrange multiplier, while the discretized 
version of Faraday's law [Eq. p3|) ] is derived as an Euler- 
Lagrange equation for the variational problem, so that 
absolute consistency with the Maxwell equations is ob- 
tained. Moreover, maximal global (integral) shielding is 
achieved through a maximal local shielding rule [Eq. [14] 
that reproduces the elementary evolution of 9tJ for a 
perfect conductor with restricted currents. 

In 3-D problems, as an advantage of the formulation 
in Eq. PH)) . one can avoid the integration of the equiv- 
alent partial differential equations and straightforwardly 
minimize the discretized integral by using a numerical 



algorithm for constrained minimization. It is this nu- 
merical minimization, instead of numerical integration of 
PDEs, which represents a very important advantage in 
the performance and power of the variational method. 
The Lagrange multiplier p (basically, the electric field) 
disappears in the direct minimization process, while the 
magnetic field can be expressed in terms of an external 
contribution /ioHo and the local sample's currents. As 
a consequence, only the unknown current components 
appear in the computation, reducing the number of un- 
known variables. Any symmetry of the problem will allow 
further simplifications and correspondingly faster numer- 
ical convergence. 

Being more specific, the integrand i(AB)^ can be 
rewritten as ^(AB) • (V x AA), and manipulated to 
get i(AA) • (V X AB) plus a divergence term, fixed 
by the external sources at a distant surface. Now, the 
integral is restricted to the samples region f2, because 
V X AB /ioAJ is only unknown within the supercon- 
ductor. In addition, the vector potential can be expressed 
as 



Mo 



AA = AAo + ^ 



A3 



47r .In k - r' 



(15) 



This transforms J- into a double integral over the body 
of the sample, i.e.: 



T= — I AAo-J„+i(r)d3r 



J„,,(r').[J„.r(r)-2J^^3,,3,,(i6) 



r - r' 



(terms independent of Jn-i-i have been omitted). 

Finally, in addition to the incorporation of the exter- 
nal sources (Aq), and the constraints for the allowed J 
region Ar one must also ensure the charge conservation 
condition by searching the minimum for the allowed set 
of current densities fulfilling V • J = as an additional 
constraint. 



III. GENERAL CRITICAL STATES: 
APPLICATIONS 

In this section, we will show that the variational state- 
ment may be used to predict the magnetic structure for 
the T-states in a three-dimensional slab geometry, i.e.: 
both in-plane and perpendicular magnetic field compo- 
nents are applied to an infinite slab and varied in a given 
fashion. A wide range of applied fields will be considered, 
and our results compared to those available in the litera- 
ture. Moreover, we will study the corrections that appear 
when a more general CT-state framework is introduced. 

First, we will give the details related to the mathe- 
matical statement of the general critical state in the slab 
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geometry. Then, the theory will be applied for estab- 
lishing the capabilities of several versions of the DCSM 

within different physical scenarios. (1 — h'l^i)^i + (1 ^ ^ '^hx^ihy^i^iipi < 



A. Infinite slab: general double critical state 

Within the infinite slab geometry, the variational for- 
mulation of the general DCSM allows an algebraic state- 
ment that is rather convenient for the eventual numer- 
ical application. To be specific, we will consider an in- 
finite slab, that is fully penetrated by a perpendicular 
uniform field H^o and then subjected to a certain pro- 
cess for the apphed parallel field, i.e.: [H^oit), Hyo{t)] 
as indicated in Fig. [5l Recalling the symmetry proper- 
ties of the electromagnetic quantities, one can describe 
the problem as a stack of current layers parallel to the 
sample's surface. Assuming that the slab occupies the 
space \z\ < d/2, it suffices to discretize the upper half, 
i.e.: < z < d/2 as symmetry (or antisymmetry) condi- 
tions may be applied. Thus, in what follows, a collection 
of N layers {zi = Si , S = d/2N) will be considered. 
Within this approximation, one has to include two com- 
ponents of J within each layer, i.e.: [Jx(zi), Jy(zi)]. No- 
tice that position independence for a given value of z 
ensures a diver genceless J. Furthermore, a sheet current 
may be introduced. Thus, the problem will be stated 
in terms of = Jx{zi)5 and -0^ = Jy{zi)S (S denoting 
the width of the layers). Now, a straightforward appli- 
cation of Ampere's law allows to express the penetrating 
magnetic field as the sums over the layers 



£.! + ijf + 2h,^,hy., e.^. < ilw ■ (19) 

(Here, the normalization h = H/iJ has been used.) 

Recall that the subindex n is introduced to indicate 
time discretization, i.e.: ipi{t+ l!^t)—ilji{t) = il'i,n+i—i^i,n- 
When this index is omitted, it will be meant that the 
equation is time independent (it is valid V n actually). 

Finally, the reader can check that a straightforward 
substitution of the squared components of the magnetic 
field entering the expression in Eq. pH)) in terms of Eq. 
([T7| leads to the following formulas for the mutual in- 
ductance coupling between layers 

Mi^. = Mf^ =1 + 2 [min {^, j}] yi^j 

A/S = A/|^^2Q+z-l) (20) 

Below, we present a number of results obtained by ap- 
plication of the previous equations for the slab. First, an 
infinite band model ( J|| ^ Jc\\ or the so-called T-states) 
will be considered. Afterward, the corrections related to 
the flux cutting limitation {Jc\\) will be studied. When 
possible, our results will be compared with available lit- 
erature. 



Hy{z,)^Hy,, = (17) 

Next, we recall that in the slab geometry Eq. (jlOp be- 
comes a discretized principle restricted to the volume of 
the slab. Following the concept introduced in the previ- 
ous section (see Eq. I16|) . the problem may be transformed 
into the minimization over the current densities 

min F = i^^,,„+iAf;5.ej,n+i-^a„M54n+i 

- ^llJi,n+l{'i - l/2)(iJj,o,n+l - H.x0^n) 
i 

+ ^6,n+l(i- l/2)(ifj,0,n+l -^f^Oai). (18) 

i 

We stress that minimization has to be performed under 
the restrictions J|| < Jc\\ and J± < Jc±, i.e.: the DCSM 
hypotheses. Specifically, one has to invoke the conditions 



1. T-states in 3-D configurations 

Here, we show the theoretical predictions for the T- 
states along the magnetization process indicated in the 
insets of Fig. [5l Starting from a fully penetrated state, 
with a field applied perpendicular to the slab surface 
(Hzo), one configures either a diamagnetic or a paramag- 
netic critical state by sweeping the applied parallel com- 
ponent Hxo (thus inducing Jy). Eventually, an increasing 
ramp in the other field component, Hyo (thus inducing 
Jx) is applied. The response of the superconductor is 
obtained as a collection of values for the sheet currents 
{Cij a-t the forward time layers n = 1,2,3,.... The 
magnetic field profiles and magnetic moments are even- 
tually obtained by numerical integration. In particular, 
the expression of the magnetic moment per unit area 

I r rd/2 
M = — / z X Jd^r = / zxJdz (21) 

JVol J -d/2 

has been used, with L representing the length of the 
sample. Recall that we have invoked the property that 
for long loops, the contribution coming from the U-turn 
at the ends, exactly equals the contribution of the long 
sides. This may be shown starting from the condition 
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FIG. 5: (color online) The magnetic moment {Mx, My) of the 
slab defined by Eq. (|21|l as a function of Hyo- Shown are 
the diamagnetic and paramagnetic cases for HxO = 1-1, and 
Hzo = 1.5 (top) and H^o — 10 (bottom). The experimental 
processes for the applied magnetic fields are shown as insets 
in the bottom panel. Units are j^x d/2 for H and ic±dV4 for 
M. Our results (lines) are compared to those by Brandt and 
Mikitik in Ref. (symbols) . In all these cases, we have taken 
Jc\\ = oo and Jc± — 1 (T-states). 



V • J = (no sources) that allows to consider the cur- 
rent density distribution as a collection of loops and en- 
sures the equality of the integrals over zJy and zJ^ (see 
Ref. Hi). 

Owing to the rich phenomenology encountered, the re- 
sults will be given separately for moderate and low per- 
pendicular fields. Recalling that H is measured in units 
of the physically relevant penetration field Jc±d/2, then 
Hzo = 0.1, 1.5, 10 will cover the range of interest. 

a. Moderate fields Fig. [5] shows our results for 
Hzo = 1.5 and H^o = 10 compared to those of Ref. 
obtained under the same conditions. A remarkable agree- 
ment is to be noticed, thus validating our theory against 
the differential equation approach of that paper. The 
same degree of coincidence was also checked for the mag- 
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FIG. 6: (color online) The magnetic moment Mx (solid lines) 
and My (dashed lines) for the T-states, as a function of Hyo- 
Shown are the diamagnetic (top) and paramagnetic (bottom) 
cases for Hzd ~ 0.1 (black), H^o ~ 1-5 (blue), and H^o ~ 10.0 
(red). Units are jc±d/2 for H and for M. 



netic field and current densities (not shown for brevity) 
as expected. In fact, the equivalence of our maximum 
projection rule and the E(J) law based analysis^ may be 
proofed as follows. The material law in that work was 
applied in two steps: (i) the transient electric field was 
chosen along the direction B x (J x B) as dictated by 
the flux flow condition, and (ii) the magnitude of E was 
found from the condition = Jci_ , through the Ampere 
and Faraday's laws and appropriate boundary conditions. 
We recall that J± = Jc± is equivalent to the selection of 
our horizontal band for Ar, that the direction of E is 
straightforwardly the same, and finally that the remain- 
ing component Jy is also coincident as it is obtained from 
the Maxwell equations, also contained in the variational 
formulation. 

b. Low fields In Fig. [6] we display the effect of 
extending the previous results to the low field region 
{Hzo = 0.1), by comparison to the values Hzo = 1.5 
and = 10. The plots indicate the following features, 
(i) in general, a saturation is reached for My{HyQ), as 
compared to the eventual linear increase of Mx{Hyo) for 
the highest values of Hyo, (ii) the higher Hzo, the sooner 
the saturation is reached, (iii) increasing Hzo rapidly di- 
minishes the slope of Mx{Hyo), (iv) in the paramagnetic 
case, a minimum is observed (more evidently for M^, and 
more visible in Fig. [5] for moderate Hzo, that is smoothed 
either for the higher or lower values of this field compo- 
nent. 

As indicated above, the underlying flux penetration 
profiles for the moderate field region were already pre- 
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FIG. 7: (color online) Profiles of the component J|| for the 
limit Jc|| — > oo (T-state) with H^o = 0.1. In all cases the per- 
pendicular current profiles satisfy J± = Jc^ ~ 1.0. The dia- 
magnetic (top) and paramagnetic (bottom) cases are shown. 
Top: inset (a) shows a zoom of J|| for the first profiles of Hyg. 
Inset (b) schematically shows the evolution of the vector J as 
function of its parallel and perpendicular components. Bot- 
tom: inset (a) shows a zoom of J|| for the first profiles of 
increasing Hyo- Inset (b) shows the magnetic moment com- 
ponents {Mx, My) per unit area as a function of Hyo. 



sented in Ref. [3 and fully coincide with our calculations. 
However, the low field region was uncovered. Here, we 
will show the peculiar behavior of the field and current 
density profiles for this regime. Thus, Fig. [7] displays the 
behavior of the projection of the current density onto the 
direction of the magnetic field (J||) under the ansatz of a 
T-state structure for H^o = 0.1. Recall that, hereafter, z 
is given in units of d/2. Then, 2 = corresponds to the 
center of the sample and 2: = 1 to the surface. 

It is apparent that the full penetration of the T-state 
requires a high field component {HyQ « 18 and HyQ w 30 
for the dia- and paramagnetic cases respectively), and a 
very high ratio J\\/Jci_ (~ 180 for the diamagnetic case 
and w 340 for the paramagnetic one). Notice that until 
these values are reached, one has J|| — 0, J± ~ 1 for the 
inner part of the sample, and a certain distribution J|| (z) 
for the outer region. We also recall a somehow complex 




FIG. 8: (color online) Magnetic field components Hx{z) (top) 
and Hy{z) (bottom) corresponding to the current density pro- 
files for the T-state limit with H^„ =0.1 (Fig. [7)). The 
diamagnetic (left) and paramagnetic (right) cases are shown. 
The insets shown a zoom of the corresponding pictures. 



structure with one or two minima in between the surface 
of the sample and the point reached by the perturba- 
tion. Interestingly, when Hyo grows, the minima become 
very flat, corresponding to a nearly constant value of J||. 
From the physical point of view, the minima basically 
represent the region where H rotates so as to accommo- 
date the penetration profile H(z) to the previous state 
of magnetization {Hj., 0, i?zo)- From the point of view of 
Faraday's law, this takes place as quickly as possible so 
as to minimize fiux variations. The obtained magnetic 
field profiles are shown in Fig. [H] Their interpretation, 
in terms of the critical current restrictions (T-states), is 
simplified by the increasing value of HyQ. Thus, a steep 
variation of Hy occurs for the inner region of the sam- 
ple, corresponding to the large values of Jy (essentially 
Jx because of the increasing in that region). On the 
other hand, displays a small slope, which relates to 
the condition Jj_ = 1 (essentially, J± ~ Jy in the inner 
region) . 

2. C'T-states in 3~D configurations 

In this section, we concentrate on the effect of con- 
sidering a flux cutting limitation (Jc||)- Magnetization 
curves, current density and field penetration profiles will 
be shown, corresponding to the same magnetic processes 
indicated in Fig. [5l but now for the rectangular DCSM 
regions with a number of values for x = Jc\\ I JcA. ■ In order 
to obtain continuity with the T-state results (recall that, 
ideally this corresponds to the limit x ~* 00) a range of 
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FIG. 9: (color online) The magnetic moment {Mx,My) of 
the slab per unit area as a function of Hyo. Shown are 
the diamagnetic (top) and paramagnetic (bottom) cases for 
HxO = 1.1 and Hzo = 1-5. The experimental conditions 
and units are those defined in Fig. [5] The T-state curves 
(Jc|| 3> 1) are shown for comparison with several rectangu- 
lar cases: J^^^ =3.0, 2.7, 2.3, 2.0, 1,7, 1.3, 1.0. The insets 
show the particular case Jc|| ~ 2.0 in the region where the 
transition T^CT is visible. 



increasing values for the parameter x will be analyzed. 
On the other hand, owing to the rich phenomenology en- 
countered, the results will be given separately for mod- 
erate and low perpendicular fields, i.e.: Hzo = 1.5 and 
H,o = 0.1. 

a. Moderate fields The main facts for H^o = 1-5 are 
shown in Figs. IMTTI First, we plot the corrections to 
Mx and My both for the dia- and para-magnetic cases, 
when the DCSM region corresponds to the aspect ratio 
values = 1.0, 1.3, 1.7.2.0, 2.3, 2.7 and 3.0 (Fig. [9|). It 
is noticeable that the limitation in J^n produces a corner 
in the magnetic moment dependencies M^.y^Hyo), which 
establishes the departure from the master curve defined 
by the T-state. The corner in and My appears at 
some characteristic field H*q that increases with x, even- 
tually disappearing within the region of interest. The 
value X* for which the corner is not observed, defines 




FIG. 10: (color online) Profiles of the parallel currents J|| for 
the T-state hypothesis "Jc|| — > oo and Jc^ = 1.0" (top) and 
for a rectangular DCSM with "J^^ = 2.0 and Jc^ = 1.0" 
(bottom). The diamagnetic (left) and paramagnetic (right) 
cases are shown. In the paramagnetic case, the profiles of Jy 
for Hyo = 1.35, 1.7, 2.0, 2.3, 2.6, 3.0 are shown as an inset, 
and correspond to the sign change in the slope of the magnetic 
moment M^. See the text for the definitions of the states 
cT'^', cT^^\ and CT. 



the range of application for the T-state limit (x^ ^ 3 in 
the conditions of Fig. [9]). On the other hand, the fine 
structure of the corner is shown in the insets of Fig. [51 
Notice that, in fact, the deviation from the master curve 
takes place in two steps, being the second one that really 
defines the corner. 

In order to allow a physical interpretation on how the 
T-states break down for the 3-D configurations studied 
in this section, in Fig. [TO] we plot the profiles of J|| (z) 
within the slab, as HyQ is increased. The upper pan- 
els show the evolution of this quantity for the T-states, 
whereas the lower panels show the process of saturation 
in which Jy reaches the value J^w both for the dia- and 
para-magnetic initial conditions. Just for convenience, 
we have introduced the following notation. cT denotes 
that J|| has reached the limit J^w only partially within 
the sample, while CT means that Jy equals Jc\\ for the 
whole range < z < d/2. For the partial penetration cT- 
states, we additionally distinguish between the so-called 
cT(i' and cT'^) phases. A s one can see m the plot, cT(i' 
means that Jy penetrates linearly from the surface until 
the limitation is reached somewhere within the sample. 
For the diamagnetic case, the profile stops at the actual 
value Jc\\ ■ However, for the paramagnetic case, the struc- 
ture is more complex. Thus, J|| penetrates linearly until 
a linear increase (toward the center) curve is reached. 
This structure is followed until the contact between both 
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lines reaches the surface. Then, the so-cahed cT*^^^ re- 
gion appears. J|| has reached J^n at the surface, and the 
whole J|| curve "pivots" around this point until the full 
CT-state is reached. We call the readers' attention that 
the initial separations of the magnetic moment from the 
T-statc master curves take place as soon as a cT^^^-state 
is obtained. On the other hand, the corners can be clearly 
assigned to the instant at which such a state disappears. 

Just for completeness, the magnetic field penetration 
profiles, corresponding to the — > oo and ^ = 1 cases 
arc shown in Fig. 1111 Notice the change in curvature 
and slope reduction in the penetration of H^, for the CT 
states. Notice also that the Hy profiles are only shown 
for the T-states, because a very similar behavior takes 
place (just differing in a small compression for the higher 
values of Hy^. 






FIG. 11: (color online) The magnetic field components Hx{z) 
(top) and Hy{z) (bottom) in the same critical states of Fig. 
1101 The diamagnetic (left) and paramagnetic (right) cases are 
shown for the T-states and the indicated rectangular regions 
of J. The profiles Hy{z) show the same behavior for oo 

(displayed) and — 2(not shown). 
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FIG. 12; (color online) Magnetic moment components 
[Mx,My) of the slab as a function of Hya. Shown are 
the diamagnetic (top) and paramagnetic (bottom) cases for 
Hxo ~ 1.1 and H^o ~ 0.1. The infinite band T-state model 
(Jc II (x) is shown for comparison with several rectangular 

cases: = 11.0, 9.0, 7.5, 6.0, 4.5, 3.0, 2.0, and 1.0. 



b. Low fields Although the general trends in the 
CT-state corrections for low H^q do not very much dif- 
fer from those at moderate field values, some distinctive 
features are worth to be mentioned for the Mx^yi^Hyo) 
curves. To start with, we recall that the corner structure 
that defines the separation of the CT curves from the 
master T-state behavior is different. Thus, as one can 
notice in Fig. [121 it is only for the higher values of 
that the separations take place abruptly. In particular, 
a smooth variation occurs for < 6 in all cases. Also 
noticeable is the change in the behavior of the initial part 
of the Mx{Hyo) curves for the paramagnetic case. Recall 
that the minimum observed for the moderate field region 
{Hzo = 1.5) has now disappeared (this can be already 
detected for the T-states). Significantly, what one can 



see as decreases is that develops a nearly flat re- 
gion at the low values of Hyo- Physically, this means that 
the initial Hx{z) profile is basically unchanged. For the 
lowest values of this can take place over a noticeable 
range of applied fields HyQ. A detail about the origin of 
this behavior can be seen in Fig. 1131 that corresponds 
to x'^ — 2. Notice the insignificant variation of as 
compared to the changes in Hy along the process. 

Also remarkable are the peculiarities of the current 
density penetration profiles for low values of Hzo- They 
can be observed in Figs. [U and [151 ttiEit reveal new phys- 
ical mechanisms, that did not appear for the moderate 
perpendicular field values. Again, the first observation is 
that the appearance of the corner in the magnetic mo- 
ment straightforwardly relates to the current density pro- 
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files. Thus, for the lower values of x (no corner present), 
the profile Jy displays a rather simple structure, basically 
jumping from to Jc\\ at some point within the sample 
(Fig. rH)) . On the contrary, for the higher values of x 
(those displaying a corner in M(HyQ) the evolution of 
the J\\{z) profiles is much more complex (Fig. [T5|) . Let 
us go into detail about these topics, part by part. 

Fig. [Til shows the profiles J\\{z) both for the dia- and 
paramagnetic cases for = 2. Recall that the evo- 
lution of the profiles with the increase of Hyo is very 
similar. The above mentioned step-like structure with 
J|| = in the inner part and J|| = J^n in the periphery 
evolves until the full penetration state J\\ = Jc\\ , V z is 
reached. On the other hand, a very interesting feature 
is to be recalled for the paramagnetic case (lower pane 
of Fig. rn)) . For the first time along the exposition of 
this paper we have met a set of conditions that produce 
an excursion of Jj_, i.e.: Jj_ = Jc± is violated during 
the process of increasing HyQ. To be specific, J± starts 
from the condition J± = Jc±, given by the initial pro- 
cess in Hxo- Then, a basically linear decrease from some 
inner point toward the surface occurs, with an eventual 
reduction to a nearly null value at some regions within 
the sample (C-states are basically provoked). Further 
increase of Hyo produces a new CT-state. This behav- 
ior is shown in a pictorial form within the insets of Fig. 
[T4l Recall that the average current density sharply tran- 
sits from a T-state {J± = Jc± , J\\ = 0) to the CT-state 
{J± = Jc± , J|| = Jc\\) for the diamagnetic case, while a 
T C — > CT evolution happens for the initial paramag- 
netic conditions. This behavior allows a physical inter- 
pretation in terms of the evolution of the magnetic field 
profiles. Thus, as stated before, the cases with small x are 
characterized by a nearly frozen profile in Hx, as shown 
in Fig. [T31 Then the structure of Hx{z) and Hy{z) is ba- 
sically a cross between two straight lines. The crossing 
point coincides with the minimum in J^{z). Recalling 



1.0 




FIG. 13: (color online) Magnetic field components Hx[z) 
(solid-lines) and Hy{z) (daslied-lines) corresponding to the 
rectangular region J^^ll = 2.0, J^^^ = 1.0 and initial paramag- 
netic conditions. The curves follow the color scale convention 
in Fig. 1101 For clarity, the Hy{z) profile corresponding to 
Hyo = 1.40 has been labeled accordingly. The analogous plot 
for the diamagnetic case strongly resembles that of Fig. [TT] 
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FIG. 14: (color online) Profiles of the parallel (J[|) and 
perpendicular (Jx) current densities for rectangular region 
"J^^ll = 2.0, Jc^ = 1.0" with H:,o = 0.1. The diamagnetic 
(top) and paramagnetic (bottom) cases are shown. Inset (a) 
schematically shows the CT structure of the full penetration 
regime in the diamagnetic CEise. The CT-C structure behavior 
of J for the paramagnetic case is shown in inset (b). 



the interpretation of the perpendicular component of the 
current density J± = dH/dz, the minima should be ex- 
pected as + Hy has a very small variation around the 
crossing point of the two families of nearly parallel lines. 

The details about the behavior of J|| for the larger val- 
ues of X are presented in Fig. 1151 that corresponds to 
the case = 7-5. Again, owing to the complexity of 
the structure, we introduce the notation cT^^^ cT^^^ and 
cT*^^', that is explained below. Let us first recall that 
the corner appears when the partial penetration regime 
cT*^^^ extinguishes and the full sample (0 < z < d/2) sat- 
isfies J± = Jc±_ and J|| = J^ii (i.e.: CT). This property 
is clearly seen in the lower panel of the figure. Thus, 
the cT(i' regime is characterized by a T region in the 
inner part of the sample (Jj^ = Jc±_ and J|| = 0), that 
abruptly becomes CT at a point that progressively pen- 
etrates toward the center (T-CT structure). At a certain 
instant, the profile becomes T-CT-T because the outer- 
most layers develop a subcritical Jy. This is called cT^^^ 
Then, the central CT band grows toward both ends. In 
first instance, the inner T region becomes CT, giving a 
global CT-T structure, that we call cT^^' . In a final step. 
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FIG. 15: (color online) Top: Profiles of J|j for the diamagnetic 
case within the rectangular DCSM with = 7.5 and H^o = 
0.1. In all cases, one gets J± = Jci_ = 1.0. Bottom: The 
corresponding magnetic moment components {Mx,My) as a 
function of Hyo are shown. The evolution from the initial 
full penetration T state to the final full penetration CT state 
takes place in three steps that are classified according to the 
structure along the sample width by: cT^^^ = T-CT, cT^^^ = 
T-CT-T, cT^^' EE CT-T and eventually CT. 



the surface T layer shrinks again to a null width and the 
full profile is a CT region. This moment establishes the 
appearance of the corner in the magnetization curves. 

B. Infinite slab: other CS models 



plifying ideas that flux depinning is only related to J± 
and that flux cutting is only related to Jy . As indicated 
in Ref. 0, in a general scenario, one should consider the 
dependencies Jc± = Jc±{J\\) and J^w = Jc\\iJ±)- In this 
work, we do not attempt a microscopic justification on 
how the DCSM hypotheses should be corrected. How- 
ever, on the basis of minimum complexity, we will analyze 
two facts: (i) the flux cutting criterion will be revised so 
as to account for the 3-D nature of the problem. In fact, 
we will show that if one considers a critical angle thresh- 
old, the cutting barrier depends both on Jy and J±. (ii) 
Also, we will investigate a smooth version of the DCSM 
in which the corners of the rectangular region have been 
rounded (see Fig. [5]). Physically, the idea behind this 
property is that the mechanisms of flux depinning and 
cutting are not fully independent after all, as one could 
expect in a continuum theory. Along this line, we recall 
that an elliptic model was introduced in Ref. [2^ that 
produces a rather good description of experimental data 
for situations in which J± and J|| occur. 

1. Critical angle gradient in 3-D configurations 

111 this part, we present some results related to the 
concept of critical angle gradient in 3-D systems. It is 
well known that, in fact, the limitation on J|| appears as 
related to the energy reduction by the cutting of neigh- 
boring flux lines when they are at an angle beyond some 
critical valuei^i^ This concept has been largely exploited 
in the 2-D slab geometry for fields applied parallel to the 
surface^ and is introduced by the local relation 



da 






dz 




H 



that establishes a critical angle gradient. Here, a stands 
for the angle between the flux lines and a given reference 
within the XY plane (i.e.: an azimuthal angle). How- 
ever, for the 3-D eases under consideration, the relative 
disorientation between flux lines may also have a polar 
angle contribution, i.e.: H does not necessarily lie within 
the XY or any other given plane. As sketched in Fig. [31 
one has to introduce the angle ip within the plane deflned 
by the pair of flux lines under consideration. After some 
vector algebra, it can be shown that, for the infinite slab 
geometry, with a 3-D magnetic field one has 

dip [ji 1 r HI ( „ TV 

(23) 

where the third component is also introduced. 

The above result is just a particular case of the relation 



As stated before, our theory will be used to investigate 
the properties related to several modifications of the con- 
ventional DCSM considered in the previous section. Such 
modifications can be justified as corrections to the sim- 



V X [Bf^ = (VB) X B 



B V X B 



[J±,l]+ [J±,2+J||] 



(24) 
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FIG. 16: (color online) The magnetic moments (solid 
lines) and My (dotted lines) of the slab as a function of Hyg 
for the critical angle gradient model [Eq. (|25|l ]. The unre- 
stricted case (kc oo) is shown for comparison with several 
cases with restricted angle gradient: k.^ =0.20, 0.30 and 0.40 
(dimensionless units are defined by Kc = Kcd/2). Shown are 
the diamagnetic (top) an paramagnetic (bottom) cases for 
= 1.1 and H^g = l.^. The insets detail the evolution of 
the angle gradient profiles for — » oo. 
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FIG. 17: (color online) The magnetic moments 
the slab per unit area as a function of Hy^^ in the diamagnetic 
(top) and paramagnetic (bottom) cases with Hx„ = 1.1 and 
i?2Q=1.5. The "infinite band" (external solid lines), rectan- 
gular (solid lines), superelliptical (dashed-lines), and elliptical 
(dotted-lines) models are shown for several values of the ratio 

>^c||/ Jcl.- 



showing that, in general, both J|| and J_l can con- 
tribute to the spatial variation of the direction B. 

Below, we display the effects of using the cutting lim- 
itation 



dz 



(25) 



instead of assuming a constant value for the parallel crit- 
ical current. Fig. [16] contains the main results. The 
calculations have been performed for the same magnetic 
processes (dia- and paramagnetic) considered in the pre- 
vious section. 

In general, (compare Figs. [9]and [T6|) one can see that 
the smaller values for the cutting threshold in whatever 
form produce the smaller magnetic moments. However, 
some important differences are to be quoted. On the 



one side, the critical angle criterion \^/\ < Kc produces 
a smooth variation, by contrast to the corner structure 
induced by the critical current one Jy < Jc\\ ■ On the 
other hand, the effect of changing the value of Kc is much 
less noticeable, especially for the diamagnetic case, in 
which the full range of physically meaningful values of Kc 
produce a negligible variation. 



We call the readers' attention that the above men- 
tioned range for Kc is established by the application of 
Eq. to the initial state of the sample. Thus, if one 
takes J|| = 0, Hzo = 1.5, H^o ~ 1.1, the squared angle 
gradient takes the value ip''^ — 0.19 and one has to use 
> 0.19 in order to be consistent with the initial critical 
state assumed. 
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2. Smooth CS models 

Here, we develop the concept of smooth double criti- 
cal state model, introduced before. Mathematically, the 
effect of rounding the corners for the rectangular DCSM 
region may be represented by a one-parameter family of 
functions with the generic form 



Jc\\ 



Jl 



< 1. 



(26) 



Such kind of curves arc known as a superellipses and cover 
the range of interest just by allowing n to take values 
over the positive integers. As the reader can easily verify, 
n = I corresponds to the standard ellipse and n > 5 is 
already a rectangle with faintly rounded corners. 

In order to illustrate the effect of smoothing the al- 
lowed region of current density components (J||, J±), be- 
low we will show the magnetization curves that arc ob- 
tained for the dia- and paramagnetic cases considered 
before. Wc compare the predictions for n — > oo, n = A 
and n = 1. For simplicity, they will be named after rect- 
angular, superelliptic and elliptic. The main results are 
plotted in Figs. [TTlfTgi 

Fig. [T7]shows the behavior of and My for an exter- 
nal perpendicular field at the moderate intensity region 
HzQ = 1.5. The first observation is that the overall effect 
of reducing the value of x = Jc\\/ Jc± is the same for the 
three models. The smaller the value of Xi the higher re- 
duction respect to the T-state (x — > oo) master curve for 
the magnetic moment components. On the other hand, 
as the particular details for each model, we recall: (i) 
the smooth models lead to smooth variations, i.e.: the 
corner is not present, (ii) the breakdown of the T-state 
behavior occurs before (at higher values of x or lower val- 
ues of Hyo) for the smoother models. Strictly speaking, 
the concept of T-state is only valid for the rectangular 
region, but it is asymptotically generated as the superel- 
lipse parameter n grows. Finally, (iii) The isotropic CS 
limit, given by the circular region n = 1 & x = 1 pro- 
duces the expected results jii collapses to zero, and 
Aly develops a one dimensional critical state behavior. 

Fig. [18] shows the comparison of and My for the 
same models considered above, but now for a low per- 
pendicular field {Hzo = 0.1). Here, one can notice: (i) 
the rectangular and superelliptical models produce very 
similar results for the diamagnetic case, both for M^ and 
for My, noticeably differing from the elliptical region pre- 
dictions, that still show a practical collapse of M^ and a 
saturation in Aly as stated before, (ii) The paramagnetic 
case involves a higher complexity. Thus, we recall that 
the already mentioned "flat" behavior of Mx for small 
values of x within the rectangular region model. This 
feature is no longer observed upon smoothing of the re- 
striction region. On the contrary, the smooth models 
involve an initial negative slope and a minimum, resem- 
bling the behavior of AI^ for the rectangular model, but 
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FIG. 18: (color online) Same as Fig. 1171 but now for HzO 
0.1. 



among the three models are also to be recalled. 

In order to provide a physical interpretation of the be- 
haviors reported in the above paragraphs, a comparative 
plot of the current density vectors for each case is given 
in Fig. 1191 For clarity, we restrict to the representation 
of the vector J at the surface of the sample {z = d/2) 
for a selected number of values of HyQ. Just at a first 
glance, one can relate the best coincidence in predicted 
magnetization to the more similar critical current density 
structures (superelliptical and rectangular regions for the 
diamagnetic case with Hzo = 0.1). Recall that, in this 
case, the rectangular region produces a CT-state struc- 
ture (J|| = Jc\\ and J± = Jci_) that is represented by a J 
vector, pinned in the corner. On the other hand, the vec- 
tor J related to the superelliptic model does not pin at 
any point, because such a singular point does not exist. 
However, it is basically oriented in the same fashion and 
this relates to the good agreement in M. We finally em- 
phasize that the cases in which strong differences occur 
for the magnetic moment, are also related to important 
changes in the behavior of J. Thus, if one considers the 
paramagnetic case at small values of Hzq and HyQ, the 
in moderate Hzo- As concerns Aly, important differences significant differences in magnetization relate to an op- 
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FIG, 19: (color online) Current density vector in the J± 
vs. J|| representation for the rectangular (red), superellip- 
tical (green), and elliptical (blue) models. The diamagnetic 
and paramagnetic cases with H^o = 1,5 and Hz = 0,1 are 
shown for several values of the ratio Jc\\/Jc± and at several 
values of the field Hyo as labeled on each arrow. Recall the 
scales on the horizontal axes that have been re-sized for visual 
purposes. 



posite behavior in J, Moreover (see left bottom panel of 
Fig, [TO)) , the rectangular model predicts a transition to- 
ward a C-state (J|| = J^w and J± « 0), while the smooth 
versions produce a tendency toward the T-state, 



C. The longitudinal problem with transport 
currents. 

In this section, we still investigate the infinite slab 
geometry, now under the assumption of a longitudinal 
transport current. Such configuration has been a long- 
standing problem, related to the design of superconduct- 
ing devices, and is still frequently focused,— and de- 
scribed in terms of the critical state regime. Here, we 
will consider the slab geometry, subjected to a uniform 
field normal to the surface (Hzo), then a transport cur- 



rent applied along the y — axis, and eventually, a mag- 
netic field {Hyo) along the same direction. This situation 
matches the second example in Ref. 0, but here, no re- 
strictions will be required for the ratios = Jcl/ Jc\\ 
and <; = Jc±d/ Hzo, that are small parameters in that 
case. Notice that the smallness of means that the 
arising critical state is approximated by the unbounded 
band region |Jj^| = Jc±,0 < |J||| < oo described be- 
fore (T-states). The smallness of ^ was meant to indi- 
cate a small deviation of the full magnetic field respect 
to the z — axis. Then, moderate values of J|| are ex- 
pected. Remarkably, these hypotheses allowed to obtain 
a set of approximate analytic formulas for the electro- 
magnetic quantities that allow to bypass the numerical 
solution of the differential equations. However, as it will 
be shown below, the range of application is narrower than 
expected. By using our numerical method, that allows 
to calculate the sample's response for any value of the 
parameters x and the range of application of such ap- 
proximation will be discussed. 

In brief, our results are not limited to the weak longi- 
tudinal current conditions. By contrast, the calculations 
are performed numerically, allowing to display the cor- 
rections needed in the general critical states. 

1. Mathematical statement 

Technically, the application of a transport current re- 
lates to the consideration of specific boundary conditions 
for the electromagnetic fields. Within our mutual induc- 
tance formulation [Eqs. (|17H20p ]. the above described 
longitudinal problem takes the following form within the 
DCSM framework. One has to minimize 

+ ^e.,„+i(*-l/2)(i?yo,n+i-i?j,o,„) (27) 

for 

hi,, it + ^2 ^ 2h,^,hy,, < Jell (28) 

and 

i^t = ^transport • (29) 

i 

This last condition, indicates that a certain transport 
current is being applied to the sample. 
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On the other hand, as related to the symmetry prop- 
erties for the transport configuration [')pi{z) = ipi{—z) 
as opposed to the antisymmetry for the case of shield- 
ing currents] , here one has to use the mutual inductance 
expressions 



Mi^. = l + 2[min{i,j}] 



Ml- = l + 2[Ar-max{2,j}] 



A/^ = 2 I - + i - 1 



Ml = 2\-+N -i 



(30) 



with the full number of layers in the discretized slab. 
As a final detail, the evaluation of the magnetic fields has 
to be made according to 



(31) 



The results obtained by application of Eqs. (|27h [3T|) 
are displayed in Figs. HD] and [211 They are described 
below. 



2. T-state solutions 
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First, we analyze the case in which the critical current 
ratio Jc±/ Jc\\ is small, i.e.: T-states are warranted. Nev- 
ertheless, here, the existence of moderate or even high 
values for the parallel component of J will be allowed. 

Fig. [20] shows the main features of our investigation. 
This includes the comparison of the penetration profiles 
for Hx and Hy obtained from our theory and from the 
analytic expressions in Ref. 0, i.e.: 



Ha; — 



Hy - 



^arcsinhfi 
cos \a 



Hyo - a 



cos 6 



2a arcsinh ( — 
a 




(32) 



Here, a has to be obtained for each value of the applied 

Hzo/ 1 



HzQ 



HyO- 



field from the condition cos9 

One can notice that the agreement is rather good for 
the higher value of H^o (200 in our dimensionless units), 
whereas remarkable differences appear for Hzo = 20 as 
HyQ increases. Our interpretation of the facts is as fol- 
lows. 

As regards the establishment of the full penetration 
profile, we have straightforwardly obtained this condi- 
tion through the step-by-step integration starting from 



FIG. 20: (color online) Profiles of the magnetic field com- 
ponents Hx{z) and Hy{z) for the longitudinal problem cor- 
responding to a transport current along the y axis of value 
/transport ~ Jc±d/2 and at several increasing values of the 
magnetic field Hyo as labeled in the curves. A slab geom- 
etry with uniform perpendicular field {H^o = 20 and then 
Hzo ~ 200) was assumed. The plot shows the comparison 
of the full range numerical solution (continuous lines) to the 
analytical approximation in Eq. (|32|l (dashed). The insets 
show the initial fiux penetration profiles for both components 
of the magnetic field. The infinite band (T-state) model has 
been assumed. 



the state HyQ = (the evolution is shown in the in- 
sets of the figure). Whereas the value 0.796 is estimated 
for the penetration field H^q within the analytical limit, 
by the straightforward method described above we get 
HyQ = 0.845. In spite of some small differences for the 
low field profiles, at moderate values {Hyo < Hzo) the 
curves always coincide. On the other hand, the failure 
of the analytical approximation for the higher values of 
HyQ is readily explained by the observation of the plot. 
Thus, increasing HyQ can compress the transport current 
toward the center of the sample (as indicated by the slope 
of Hx{z)). For the case of H^q = 20, one gets Jy.max ~ 5 



when Hy w 100 and Jy 



50 when Hy « 1000, then 
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a considerable value of J\\ is obtained. This leads to a 
not so good approximation from the analytic condition 
in the approximation of Ref. 0, which one is only valid 
for small values of this quantity. However, when compar- 
ison is made for Hzo = 200, one gets Jy^max ~ 1 when 
Hij « 100 and Jy.max ~ 5 when Hy w 1000. Then, a 
much better performance is obtained for the analytical 
limit even for very high applied fields HyQ. 



3. CT-state solutions 

This part will be devoted to unveil the features of lon- 
gitudinal transport problems under general critical state 
conditions. To be specific, we will compare the transport 
current profiles obtained either by assuming a T-state 
(flux cutting may be neglected) or a general critical state 
in which flux pinning and transport are at the same level, 
i.e.: X = Jc\\/Jc± = 1- Pictorially, (see Fig. [^, we 
solve the CS problem either for a horizontal band or for 
a smoothed square region. All the results shown corre- 
spond to the case H^q = 20. 

First, let us recall that the initial transport profile 
(Jy(z) = for z < 0.5 and Jy{z) = 1 for z > 0.5) is 
spread out by the action of HyQ. Thus, as illustrated 
in the upper panel of Fig. [211 the free current core is 
reduced stcp-by-step until the eventual full penetration 
occurs for an applied field HyQ = 0.845. We remark that 
no significant difference is observed when this plot is gen- 
erated, either for the T or CT critical state models. On 
the other hand, the full penetration regime displays clear 
differences, at least for the high field region. Thus, for 
the case under consideration (iJzo = 20) one can no- 
tice that the transport current profiles are very similar 
for HyQ < 20, but display even qualitative differences for 
Hyo > 20. This property relates to the appearance of the 
limitation for J||. In fact, the value J|| = 0.9 Jc\\ = 0.9Jc_l 
is obtained at z = when HyQ equals 20. Subsequently, 
J|| (which one can basically identify with the transport 
current Jy for large values of Hy) increases more and 
more at the central region of the sample for the T-state 
model. For the CT-state, the situation at high fields is 
rather different. Initially, Jy reaches the limit Jf,|| (= Jcj_ 
in the case under study) at the center, and rapidly de- 
creases toward the value JcA_/2, that is roughly main- 
tained for all z > 0. 

The behavior of J is obviously inherited by the flux 
profiles. One can see it in the plot of H^ (see the insets 
of Fig. [^T|) . For the T-states, H^ saturates at the periph- 
ery, where a practically fiat profile is reached. Then, one 
has ^ in that region, and the full transport profile is 
shifted toward the center of the sample. At the same time, 
the shielding part of the critical current « J± is max- 
imum within the region of negligible transport and goes 
to zero where transport predominates. Such a behavior 
(compression of transport by penetrating shielding cur- 
rents) is straightforwardly deduced from the slopes of the 
magnetic field profiles in Fig. [201 and has been already 



suggested as a possible explanation the magnetic field 
dependence of the transport critical current in longitu- 
dinal geometryi^ Direct measurements of the transport 
current density profile;^ have also been used to conclude 
that a longitudinal field compresses the transport current 
toward the center of the sample. On the other side, for 
the CT-state case, in which the limitation on the parallel 
current density is active, H^ displays a turn back until 
a nearly linear penetration is reached, basically charac- 
terized by the slope Jc_l/2. Thus, the transport profile 
eventually stretches instead of concentrating toward the 
center. 



IV. CONCLUSIONS 

In this article, we have shown that the critical state the- 
ory for the magnetic response of typc-II superconductors, 
may be built in a very general framework. The basic con- 
cepts underlying the phcnomcnological approach, issued 
by C. P. Bean in the early sixties have been identified as 

1. The CS theory bears a Magneto Quasi Stationary 
approximation for the Maxwell equations in which 
E and p arc second order quantities. This means 
that the magnetic flux dynamics is described by 

[ AB = V X (EAt) (implicit) J 

i.e.: the inductive part of E may be introduced 
through Faraday's law, whereas the role of electro- 
static quantities is irrelevant. E may be modified 
by a gradient (E E + V0) with no effect on the 
magnetic response. 

2. The law that characterizes the conducting behavior 
of the material may be written in the form 

[ li Ej^Q ^ max J • E J 

In physical terms, the material "reacts" with a max- 
imal shielding rule when electric fields are induced. A 
perfect conducting behavior characterizes the magneto- 
static equilibrium when external variations cease. In all 
cases, J is constrained within some region Ar. 

We stress that the importance of E is sometimes veiled 
by the actual application of the above rules, as it plays 
an implicit role usually. 

On the other hand, the above representation may be 
understood as the macroscopic counterpart of the under- 
lying vortex physics. Thus, the physical barriers for flux 
depinning and cutting are represented by the condition 
J G Ar, i.e.: the current density is confined within some 
region Ar (J < Jc in 1-D). The evolution from one mag- 
nctostatic configuration to another occurs through the 
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FIG. 21: (color online) Profiles of the transport current density Jy{z) for the longitudinal transport problem in the slab for 
the same conditions of Fig. [20] (here H^o = 20). The upper panel shows the initial partial penetration process induced by 
application of increasing values of Hyo as labeled in the successive curves. The transport current profile fully penetrates at 
Hyu ~ 0.845. The lower panels show the evolution of Jy{z) from the fully penetrated state Hyo > 1.0 both for the T-state 
model (left) and for a superelliptical model with n = 4 and x = 1 (right). The induced magnetic field profiles Hx are shown as 
insets. 



local violation of this condition, i.e.: J ^ Ar (J > Jc 
in 1-D problems). However, owing to the high dissipa- 
tion, an almost instantaneous response may be assumed, 
represented by a maximum shielding rule in the form 
max J • E Ijg^^ (J = ±Jc in 1-D). 

The general CS theory exposed above may be solved 
in different forms. In our work, we emphasize the per- 
formance of variational methods for solving the problem. 
In particular, the mutual inductance representation with 
J(r) as the unknown, offers two important advantages: 
(i) intricate boundary conditions and infinite domains are 
avoided, (ii) the transparency of the numerical statement 
and its performance (stability) are outlined. Thus, the 
quantities of interest (flux penetration profiles and mag- 
netic moment) are obtained by integration and additional 
smoothing is ensured. To be specific, upon discretization, 
the CS problem bears the algebraic expression 



(I \ 

min F = - X^ij -^4,n+lMy ^j,n+l 

with {/i,n+i} the set of unknown current values at the 
specific circuits for the problem of interest, M^- their mu- 
tual inductance coupling coefficients, and AHi the ap- 
plied magnetic field increment. Corresponding to the CS 
rule J 6 Ar, each value must be constrained. Also, 
we have found that a number of constraints related to 
physically meaningful CS models may be expressed in 
the algebraic form 
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with /o some constant representing the physical thresh- 
old, Fa an algebraic function representing the physical 
model, and C°j a coupling matrix, that also depends on 
the model. For example, the isotropic model corresponds 
to F{x) = X ; Cij = Sij ; /q = j^. On the other hand, 
the double critical state model is given by 

Fi{x) = X ; -^2(2;) = X 
foi = jI± ; /02 = Jc|| 

and the coupling coefficients C^^'^ that project the local 
current density onto the local magnetic field or its normal 
plane, are obtained from Eqs. (fT7|) and (fT9|) . 

General critical state problems have been solved for a 
number of examples, within the infinite slab geometry. 
All of them share a 3-D configuration for the magnetic 
field, i.e.: H = (Hj., Hy, H^) under various magnetic 
processes and models for the critical current restriction. 
Thus, we have considered several physical scenarios clas- 
sified by the ansatz for the flux depinning and cutting 
processes (basically affecting the critical current thresh- 
olds Jc±_ and Jc||) and their relative importance (given 
by X = Jc\\/Jc±)- In summary, the following cases have 
been analyzed: 

1. T-statc solutions, in which the approximation 
X ^ 1 produces the result J± = Jc± and J|| 
may be arbitrarily high. Our predictions show an 
excellent agreement with previous results in the lit- 
erature, and extend the theory to the full range of 
applied magnetic fields. 

2. CT-state solutions in which x Jc] 1 for several cases 
within the rectangular region given by J± < Jc±_ 
and J|| < Jc\\ are predicted by the theory. Out- 
standingly, the appearance of the flux cutting lim- 
itation takes place as a sudden corner in the mag- 
netic moment curves in many cases. The corner 
establishes a criterion for the range of application 
of T-state models. 

3. The critical angle (between vortices) criterion that 
establishes the limitation on J|| has been modifled 
for 3-D problems. It is shown that, in general, the 
concept may involve both J|| and J± as one can see 
in Eqs. ipfeS]) . 

4. The possible coupling between the flux depinning 
and cutting limitations has been studied through 
the solution of smoothed DCSM cases. In particu- 
lar, we consider the effect of rounding the corners 
of the rectangular region J± < Jc± and J|| < J^n, 
by the superelliptic region criterion {J\\/ Jcw)"^"" + 
{J^/Jc^)"^^^ < 1 with 1 < 77, < 00. It is shown 
that, under specific conditions (paramagnetic ini- 
tial state and low perpendicular fields), important 



differences in the predictions of the magnetic mo- 
ment behavior are to be expected. The differences 
in M have been related to the behavior of the criti- 
cal current vector Jc around the corner of the rect- 
angular region. 

5. The longitudinal transport problem, i.e.: a mag- 
netic field is applied parallel to the transport cur- 
rent, has been studied for several 3-D configura- 
tions. It is shown that the transport current is es- 
sentially compressed toward the center of the sam- 
ple by the effect of shielding currents when no lim- 
itation on J|| is active (T-states). However, in- 
creasing the parallel field when the constraint Jf,|| 
is reached, produces a flattening on the transport 
current density that becomes nearly uniform across 
the sample. 

We emphasize that the scope of our theory is rather 
beyond the actual examples treated in this article. On 
the one side, we have shown that the CS concept allows 
arbitrariness in the presence of electrostatic charge and 
potential, and one could simply upgrade the models by 
the rule E ^ E 4- V(f> if necessary. For instance, a scalar 
function (j) may be introduced if the direction of E has 
to be modified respect to the maximum shielding rule in 
the MQS hmit. 

On the other side, the extension of the theory to ar- 
bitrary sample geometries is intrinsically allowed by the 
mutual inductance representation. This article has laid 
necessary groundwork for attacking general critical state 
problems in 3-D geometry. Experimental studies that 
could reproduce the situations considered in the differ- 
ent examples, as a means of testing the double critical 
state model predictions are suggested. 

From the theoretical point of view, a relevant technical 
issue to be considered is that the divergenceless charac- 
ter of the current density is not always ensured. In this 
work, the problem's symmetry has allowed to identify the 
elementary current circuits that fulfill such condition (in- 
finite horizontal layers), and the corresponding geomet- 
rical problem of finding their coupling matrix coefficients 
has been solved. In general, this is not a trivial issue 
and one will have to incorporate the additional restric- 
tion V • J = or use a representation with appropriate 
basis functions for Ji^ 
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